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SUMMARY 


An analysis of the relative frequency of occurrence of any given 
slip-line angle in a plastically deformed polycrystal composed of face- 
centered cubic crystals is presented for the case of simple tension. 

The results are compared with those obtained for a polycrystal com- 
posed of crystals which have but a single mode of slip and with experi- 
mental results. The comparisons show that the differences between the 
results obtained by the two theories become greater as the stress is 
increased. The comparison of the face-centered cubic theory with 
experiment is somewhat better than that of the single -slip-mode theory, 
but the errors are appreciable. 


INTRODUCTION 


The frequency distribution of the angular orientation of slip 
lines that are observed within separate grains on the surface of a 
plastically deformed polycrystal depends upon the detailed mechanism 
of plastic deformation. In reference 1 an attempt was made to assess 
quantitatively the assumptions on which the slip theory of plasticity 
(ref. 2) is based by investigating the implications of this theory 
concerning this frequency distribution. This assessment was made by 
comparing an experimental distribution with theoretical distributions 
calculated on the basis of the same model as that used in formulating 
the stress-strain laws of the slip theory. Although good agreement 
was obtained with regard to the shape of the distributions, the com- 
parison between the experimental maximum slip angle and that predicted 
by theory was poor. One of the possible reasons for this poor com- 
parison, as reported in reference 1, was the fact that the theory was 
based on a polycrystalline aggregate of grains which possess only one 
mode of slip; whereas aluminum, the metal used in the experiment, is 
made up of face-centered cubic crystals which have 12 modes of slip. 

In order to investigate the quantitative effect of the multimode 
property of face-centered cubic crystals on the slip-angle distribution 
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an investigation was performed wherein this distribution was derived 
on the basis of the same assumptions as those used in reference 1 
except that the grains were assumed to be face-centered cubic crystals 
instead of the single-slip-mode type. The theoretical derivations 
and the results of the investigation are presented herein. 



SYMBOLS 

N(e) 

cumulative probability of slipped grain having slip angle 
less than 9 

s(e) 

probability density of slipped grains having slip angle of 6 

G(0,\) 

indie ial cumulative probability 

K(e,x) 

indicial probability density 

w(\) 

weighting function (see eq. (7)) 

e 

slip angle 

®max 

maximum slip angle 


( 3 , 03,5 coordinates specifying orientation of slip-plane-direction 



combination with respect to specimen axis and viewing 
plane (see fig. 3) 

x,y,z 

coordinate axes (see fig. 3) 

a 

tensile stress 

°L 

lowest value of tensile stress to cause slip 

R 

stress ratio, o/a l 

T 

resolved shear stress 

T L 

limit shear stress 


THE FACE-CENTERED CUBIC CRYSTAL 


The structure and distortion of the face-centered cubic crystal 
have been fully discussed in the literature (see, for example, refs. 3 
to 5). For convenience, however, some of the characteristics of single 
crystals are given here. 
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Figure 1 indicates the basic slip systems of the crystal. In 
terms of the Miller indices (which are merely sets of direction num- 
bers referred to the edges of the cube), the normals to the planes of 

slip are of the form {lli^ (i.e. (ill), (ill), (ill), etc.) and the 
slip directions in the plane are of the form <110> (i.e. [JLICJ] ^ [no], 

CloiJ, etc.). The octahedral planes therefore are the slip planes and 
the face diagonals are the slip directions. The regular tetrahedron 
shown inscribed in the cube contains the slip planes as faces and the 
slip directions as edges. Four slip planes and three slip directions 
in each plane produce a total of 12 possible plane-direction combina- 
tions or slip systems. 

If a tensile stress a is applied to a single crystal specimen 
along the specimen axis the resulting resolved shear stress on any 
particular slip system is given by 


r = a cos X cos 5 (l) 

where X and S are the angles between the specimen axis and, 
respectively, the normal to the slip plane and the slip direction. 

This shear stress is, in general, different for each slip system in 
the crystal. The particular slip system for which t has a maximum 
magnitude is, of course, dependent on the orientation of the crystal 
and can be found by comparing the shears on the various slip systems. 

The stereographic projection of the crystal shown in figure 2 summarizes 
the information obtained in reference 5 by such a comparison. The 
crystal is presumed to be fixed and the orientation of the crystal with 
respect to the loading direction is specified by the position of the 
specimen axis on the stereographic projection. The projection shows 
the orientation of the slip planes (of the form "{ill} > indicated by 
the symbol A) and the slip directions (of the form <110>, indicated 
by the symbol o) with respect to the orientation of the cube edges 
(of the form <L00>, indicated by the symbol □). The great-circle 
arcs connecting the various directions divide the hemisphere into 
2k equal triangular regions. The set of indices within each region 
expresses the particular slip plane (upper indices) and slip direction 
(lower indices) most highly loaded if the specimen axis falls within 
that region. Since the resolved shear stress required to cause slip 
along any given slip system of a virgin crystal is independent of the 
orientation of the slip system, the sets of indices also denote the 
particular slip system which will initially undergo slip. (An excep- 
tion to this statement is noted in ref. 6 wherein the required resolved 
shear stress is found to be different for different slip systems in a 
thin, plate-like crystal. ) As the crystal is loaded further in tension 
the same slip system will continue to undergo slip provided the crystal- 
lographic orientation remains essentially unchanged with respect to the 
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loading direction; the other systems will remain dormant because, as 
has been shown experimentally, the dormant slip systems are strain- 
hardened to at least the same extent as the operative ones. (See 
ref. 3 , p. 304. ) 


THEORY 


In the derivation a polycrystalline specimen is considered to be 
plastically strained in tension. The grains, which are so oriented 
that the resolved shear stress on the most highly loaded slip-plane— 
direction combination is sufficiently great, exhibit slip along the 
slip plane which is evidenced in the form of slip lines on the surface 
of the specimen. If the specimen axis is vertical, the slip lines are 
inclined to the horizontal at various angles known as slip angles. As 
in reference 1 the objective is the calculation of the relative frequency 
of occurrence of slipped grains with various slip angles. 

The desired result can be expressed, in one way, as a cumulative 
probability, that is, as the relative number of slipped grains having 
slip angles less than a given value. This cumulative probability N(0) 
can be found conceptually by counting the slipped grains with slip 
angles less than 0 and dividing by the total number of slipped grains. 
Another way to express the desired result is as a frequency distribu- 
tion or probability density, that is, as the relative density of slipped 
grains with a given slip angle. The probability density S (0) is 
found to be equal to the derivative of N(0) with respect to 0. Both 
of these quantities are derived herein. 

The following assumptions, which with the exception of assumption 2 
are common to those used in reference 1, are used in the derivation: 

(1) The crystallographic orientation of the grains in the specimen 
is random. 

(2) The specimen is composed of face-centered cubic crystals. 

(3) The microscopic stress state in each grain in the same as the 
macroscopic stress on the specimen as a whole. 

(4) A grain slips when the resolved shear stress in the slip direc- 
tion in the slip plane is greater than a certain limiting value, herein 
called the limit shear stress, which is the same for all grains. 

The model thus consists of a specimen composed of a very large 
number of randomly oriented grains of identical crystallographic pro- 
perties subjected to a common tensile stress in the direction of the 
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specimen axis. In the determination of the slip-angle distribution, the 
following questions must he answered for each particular grain: 

(1) Which is the most highly loaded slip system? 

(2) Does this system undergo slip? In other words, is the resolved 
shear stress greater than the limit shear stress? 

(3) If slip occurs, what is the slip angle produced? 

With these questions answered for each grain the relative frequency of 
occurrence of slipped grains with the various slip angles can be deter- 
mined. Now, the concept of having many randomly oriented grains can be 
replaced by the concept of having a single grain which is allowed to 
assume randomly all orientations; each different orientation represents 
a different grain. Furthermore, since the 12 slip systems are crystal- 
lographically equivalent, attention can be directed to only tho^e 
orientations for which a particular slip system (say the (lll/Oll) sys- 
tem) would be most highly loaded. If the orientation is such that 
another slip system is most highly loaded, the basic tetrahedron of 
figure 1 can be rotated to bring the (lll/Oll) system into coincidence 
with this most highly loaded system and this new orientation will be 
crystallographically the same as the old. By means of this artificial 
restriction of the orientation the first question is automatically 
answered and only the second and third questions need to be answered. 

In order to answer these questions, the orientation of each crystal 
with respect to both the specimen axis and the viewing plane must be 
specified. One method of specification was used in reference 1 wherein 
the specimen axis and viewing plane were fixed and the orientation of 
the grain was allowed to vary. A better method for the purposes of 
this analysis would be to hold the grain fixed and allow the orienta- 
tion of the specimen axis and viewing plane to vary. The quantities 
necessary to specify these orientations are shown in figure 3 wherein 

the (ill) plane is the xy-plane and the [oil] direction is the x-axis. 

The orientation of the specimen axis is given by the spherical coordi- 
nates A,, the angle between the specimen axis and the normal to the 
slip plane, and 3, the angle between the slip direction and plane OAC 
which passes through the specimen axis and is normal to the slip plane. 
The orientation of the viewing plane (which is, of course, parallel to 
the specimen axis) is given by a), the complement of the dihedral angle 
between the viewing plane OAB and the aforementioned plane OAC. The 
slip line produced by the intersection of the slip and viewing planes 
and the slip angle 0, the angle measured in the viewing plane between 
the slip line and the normal to the specimen axis, are also shown in 
figure 3. These quantities defining the orientation are exactly 
equivalent to those used in reference 1. 
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The value of co is in no way involved in answering the question 
of whether a grain will slip; variation of oi merely changes the slip 
angle produced and only the values of X and 3 need to be known to 
answer this question. If the stereographic projection shown in figure 2 
is rotated to bring the normal to the (ill) plane into the center, the 
projection in figure 4 is produced. In figure 4 only those parts of 
the projection which will be useful are shown. As was previously 
remarked, the only orientations of the_grain that are considered are 
those for which the slip system (lll/Oll) is most highly loaded. The 
triangular region specifying orientations of the specimen axis (or 
values of X. and 3 as indicated in fig. 4) for which this slip sys- 
tem is most highly loaded is therefore shown. Thus, the restricting 
condition is merely that only values of X and 3 within this triangle 
are to be considered. This limitation entirely takes into account the 
multimode property of the face-centered cubic crystal and constitutes 
the significant difference between multimode and single-mode crystals; 
there would be no restriction of values of X and 3 for a- single- 
mode crystal. 

The question of whether a grain will slip can now be answered. 

The resolved shear stress on the (lll/Oll) slip system due to a tensile 
stress a along the specimen axis is given by 


t = a cos X sin X cos 3 (2) 

The locus of positions of the specimen axis for which t = the limit 
shear stress, is given by 


t L = o cos X sin X cos 3 


or 


= sin 2X cos 


3 


(3) 


where R 


is the stress ratio and is equal to 


a 

a L 


in which 


a L = ^L- 


For a fixed value of R, equation (3) describes a closed curve such 
as that shown in figure 4 which separates the sphere into two regions, 
the one inside being where t > tl and the one outside being where 

t < tl. As R is increased the curve encompasses an ever widening region 
starting with a point at X = -j^, 3 = 0 for R = 1 and increasing to 
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the whole region 0 < X. < J, - | < 0 < | f or R = oo. Such curves and 

regions have been constructed by Von Goler and Sachs (see, for example, 
ref. b, p. 38 ) for a different orientation of the basic triangle than 
the one contained herein. 

For a grain to slip at a given stress ratio the grain must be 
oriented so that the specimen axis falls in this region. In addition, 
however, a restrictive condition is that the specimen axis must fall 
within the spherical triangle. The specimen axis thus must be included 
both in this region and the spherical triangle. For small values of R 
the curve for r = t l falls entirely within the triangle and the results 

should be the same as those obtained in reference 1 because the restric- 
tive condition is not involved. For large values of R the curve lies 
entirely without the triangle and the exact value of the stress ratio 
is unimportant; the value of R is high enough to cause slip in all 
the grains and the slip-angle distribution is frozen. At intermediate 
values of R both the triangular boundary and the value of the stress 
ratio are important. 

Attention up to this point in the discussion has been directed 
mainly to considering the variables X and 0 . Variation of cjd, as 
was pointed out, only changes the slip angle produced. The effect 
of ca is taken into account in the following manner. Define the 
indie ial cumulative probability G(X.,0) as the probability of a slipped 
grain, the orientation of which with respect to the specimen axis is 
given by X and 0, exhibiting a slip angle less than 9. (This indicial 
cumulative probability is independent of 0 , as will be seen later.) 

The equation for N(9) can then be obtained by superposition: 


where the quantity sin X, dX, d0 is an elemental area on the sphere 
shown in figure !<•, and the integrations in both the numerator and the 
denominator are carried out over the area A which is included both 
within the spherical triangle and the curve t = t^. The numerator 

thus is effectively the sum of the indicial cumulative probabilities 
for each position of the specimen axis (X,0) over all positions which 
produce slip, and the denominator is effectively the sum of all posi- 
tions (X.,0) which produce slip. 


N(0) = 


JI G(X ' 0) sin X dX dp 



oo 
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In order to derive the expression for G(X,0), the relationship 
between the orientation of the slip system (given by (X.,3,03, )) and 9 
must be known. This relationship can be obtained from figure 3 hy 
solving the irregular tetrahedron OABC and is 


sin cn 


tan 9 
tan X 


(5) 


Note that the value of 3 does not enter into this relationship since 
the slip angle produced by slip on a given system depends only on the 
orientation of the slip plane and not on the direction of slip within 
the plane. Because of the symmetry of this equation only values of 03 
between zero and ir/2 need be considered in the calculation of G(X,0). 
The quantity G(\,0) is by definition the proportion of all values 
of cn in this range for which the slip angle is less than a given 
value of 9. When X < 9, the slip angle produced will be less than 9 
for all values of oo; when X > 9 , the slip angle produced will be less 
than 9 for all values of 03 between 0 and sin-l(tan 9 /taxi X). Thus, 


G(X,0) = 1 


p . -1 tan 9 

— sin 

jt tan X 


{X < 9) 
(X > 9) 


( 6 ) 


Note that the indie ial cumulative probability is independent of 3 as 
a consequence of the fact that 3 does not appear in equation (5). Since 
the limits on the 3 -integrations in equation (4) are functions only of X , 
the expression for N(0) can be written 


N(9) 



W(x)G(X,0)sin X dX 



W(X)sin X dX 


( 7 ) 


In this expression, the term W(x) is merely the difference between the 
values of 3 at the end points of the small circle with radius X which 
is included in both the triangular area and the region inside the 
curve r = t^. The value of W(X) is zero for values of X for which 

there are no positions of the specimen axis within both of these regions. 
The expression for the probability density S(0) can be obtained by 
taking the derivative of N(0) with respect to 9 
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s(e) 



W(X.)K(X,,0)sin X. dX, 



W(x) sin \ dX 


( 8 ) 


where 


K(x,e) = ae G{X,e) 


= o 


sec 


20 


\J tan2\-tan20 


(X < e) 
(x > e ) 


S (9) 


In equations (7) and (8), the function W(X,) can be interpreted as 
being a weighting function used in the superposition of the indicial 
cumulative probability G(X,0) and the indicial probability density 
K(X,0). In order to determine the end-point values of P for the 
purpose of calculating W(x), vise must be made of equation (3) as well 
as the equations of the great circles bounding the triangular area; 
these equations are obtained in the appendix and the results are summa- 
rized here. (Each great circle is denoted by the points through which 
it passes. ) 

For great circle (001,101), 


3 = £ - cos-1 H21A (10a) 

6 \f2 

For great circle (001,111), 


3 


COS"l (/2 COt X) - y 
v 6 


(10b) 
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For great circle (101,111), 




(10c) 


As can be expected, for a fixed value of R, the variable 


X must 

be separated into various ranges in the determination of W(x). In 
addition, the stress ratio itself must be separated into various ranges. 
A discussion of each range of R and illustrative sketches follow: 


(l) Range I: 1 < R < 


JL 


i + 


/2 



In this range the curve for t = t l falls entirely within the triangle; 

the upper limit of this range is the value of R for which the 
curve r = t l first touches the great circle (001,101). The values 

of W(X) are 

W(k) =0 (0 < X < X]_) 


= 2 cos 


= 0 


-1 


R sin 2X 


(A.^ K. X < X 2 ) 
(X-2 < X ^ 


> (ID 
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where the minimum and maximum values of \ on the curve T = t l are 




— sin 
2 





Tt 

2 


1 

2 


sin 


-1 


1 

R 


(2) Range II: 


\fe 

< 

1 + \J2 


R < 


5 

2V^ 



In this range and in subsequent ranges both the curve and the triangle 
are involved. The upper limit of this range is the value of R for 
which the curve intersects the great circle (001,101) at P = 0. The 
values of W(\) are 
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w(x) = o 


= 2 


COS 


-1 


1 

R sin 2X 


= cos 


-1 


cot X 


R sin 2X 


+ cos 


i2 


Jt 

6 


= 2 cos' 


R sin 2X 


= 0 


(0 < X < X]_) 
(Xi < X < A]_) 

(A^ < X < A2) 

(a 2 < x < x 2 ) 

^x 2 < x < f ) 




( 12 ) 


where the values of X at the intersections between the curve for 
t = t l and the great circle (001,101) are 


A 


1,2 


1 

2 


cos 


K 

R 


2 ± 



3 


_ 6 _ 

R 2 


5 6 

—7c < R < p 

2/6 3 + \/6 


(3) Range III: 
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The upper limit of this range is the value of R for which the curve 
first touches the great circle (001,111). The values of W(X) are 


(0 < X < Ai) 


(A 1 < X < A 2 ) 

> (13) 

(A 2 < X < X2) 


(X 2 < X < |)^ 



w(x) = 0 


_1 1 _1 cot X 

cos ^ + cos 


R sin 2A, 


/ 2 - 


2 cos"-*- i 

R sin 2X 


= 0 


(4) Range IV: < R < 11 


3 + V /6 4 \[6 
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The maximum value of R in this range is that for_which the intersec- 
tion between the curve and the great circle (001,111) occurs at 3=0. 
The values of W(X) are 


W(\) = 0 


(0 < X. < A x ) 


i 1 -l cot X 

cos -l + cos — 


R sin 2X 


/? 6 


(a^ < X < A2) 


= 2 cos 


-1 1 
R sin 2X 


cos - cos~^(\/2 cot x) + 

R sin 2X 


.-1 L 


(A2 < X < B]_) 


(B x < X < B 2 ) 


> tut) 


= 2 cos' 


R sin 2X, 


(Bg < X < X 2 ) 


= 0 


C 


Xo < A. < 


where the values of X at the intersections between 
t = t-^ and the great circle (001,111) are 


the curve for 
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( 5 ) Range V: 


ii- < R < — 

k\/6 /3 



The upper limit of this range is given by the value of R for which 
the curve first touches the great circle (101,111). The values of W(X) 
are 


w(x) = 0 


= COS" 


= 2 cos 


1 -1 cot X 

+ cos”- 1 - 

R sin 2X ^2 


-1 L 


n 

6 


R sin 2X 


= cos -1 cos 1 (|/2 cot \) + 4 

R sin 2X ° 


= 0 


(0 < X < Ap) 
(Ap < X < a 2 ) 


(A 2 < X < Bi) S ( 15 ) 


(Bp < X < B 2 ) 


(b 2 <x< I) 


1 6 
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(6) Range VI: 





The maximum value of R in this range is that for which the curve 
passes (simultaneously) through both the corners (001) and (101). The 
values of W(X) are 


W(X) = 0 


= cos 


-1 


= cos 


= cos 


R sin 2X 
-1 cot X 

-1 1 _ 

R sin 2X 


+ cos 


_]_ cot X 


+ cos 


.1 cot X 

\[2 


= 2 cos -1 


R sin 2X 


= cos 


-1 


R sin 2X 


cos" 1 (^2 cot x) + 


(0 < X < Ap) 

(A x < X < Ci) 

(C x < X < C 2 ) 

(C 2 < X < Ag) y, ( l6) 

(Ag < X < Bi) 

(Bi < X < Bg) 
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where the intersections between the curve for t = and the great 
circle (101,111) are 


Cl 


1 

2 


sin 


2 

73R 


C 2 


it 

2 


1 

2 


sin 


2 


(7) Range VII: 




The upper limit of this range_ is the value of R for which the curve 
passes through the corner (ill). The values of W(\) are 
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W(X) = 0 


= COS 


cot X 




(o< \ < cos- 1 

(— 1 /|<x < 008-1 ^/f) 


It -1 

cos 

3 


(1/2 cot x) 


(c-- 1 ^| < x < c 2^ 


/> (IT) 


cos -1 - COS -1 ^1/2 cot x) + 

R sin 2X 6 


= 0 


(8) Range VIII: 




(C 2 < x < b 2 ) 


(b 2 < X < I) 
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In this last range the curve lies outside the triangle and the exact 
value of R is unimportant; slip occurs in all grains and the slip- 
angle distribution is frozen. The values of W(x) are 


W(X) = 0 


= cos 


-1 CQ t ^ 


= — - cos 


72 

~^(^2 cot 



= 0 



The theoretical derivation is thus completed. By using equations (ll) 
to (18) in conjunction with equations (6) and ( 9 ), the cumulative prob- 
ability N (0) (eq. (7)) and the probability density S(0) (eq. (8)) 
can be evaluated. 


RESULTS AND COMPARISONS 


In order to show the difference between the results_obtained 
herein and those derived in reference 1 the curves for N(0) and S(0) 
have been computed by both theories for three values of R and compared 
in figures 5 to 7. In the calculation of these curves, the integrations 
necessary for computing S (0) were carried out numerically. For 
isolated cases, however, the integrals could be expressed in the form 
of elliptic functions and these exact evaluations were used as a check. 
The curves for N(0) were found by numerically integrating the curves 
for 3(0), S (0 ) being the derivative of N(0) with respect to 0. 

As can be seen from figure 5 for which R = 1.015, there is no plottable 
difference between the results for this value of stress ratio. When 
the stress ratio is increased to 1 . 15 6 , the effect, as shown in figure 6, 
increases appreciably. Even more effect is apparent in figure 7, for 
which the stress ratio is 1.837* It should be noted that the value of 
the stress ratio R = 1.837 is the one for which the slip-angle dis- 
tribution derived in the present paper becomes frozen. For higher 
values of the stress ratio the curves for the face-centered cubic dis- 
tributions would be the same as those shown in figure 7> whereas the 
curves derived on the basis of single-mode crystals would continue to 
change . 


20 


NACA TN 2777 


The manner in which the differences between the two theories exhibit 
themselves can be seen from figures 5 to 7- The values of S(0) for 
the single-mode crystal remain almost constant for low values of 9 as 
the stress ratio is changed; whereas the corresponding values for the 
face-centered cubic crystal decrease as R is increased (until R = 1.837> 
of course). The value of slip angle where the probability density is 
maximum continually decreases as the stress ratio is increased for the 
single-mode crystals; whereas the value for which the maximum is obtained 
remains at about h-5° for the face r centered cubic crystal. The value of 
the maximum slip angle, denoted herein by 0max> Tor the face-centered 
cubic crystal is always less than or equal to ®max Tor the single-mode 
crystal. 

More detailed consideration can be given to this last point. It 
was found in reference 1 that ®max i s > Tor single-mode crystals, 
equal to \ 2 . The value of f^ax as obtained herein is clearly given 
by the maximum value of X, for which W(k) is not zero. (See eqs. (8) 
and (9).) Thus, for face-centered cubic crystals. 


®max — ^2 



= B 2 



<R<| 



= cos 


-1 1 




) 


These results for 0 max are plotted against the stress ratio R in 
figure 8. For comparison the corresponding results for the single- 
mode theory are also shown. 

In reference 1 an experimentally obtained slip-angle distribution 
was shown. This distribution was taken from a photomicrograph of a 
polished 2S-0 aluminum alloy specimen at 0.022 strain. At this strain, 
the stress ratio was found to be approximately 2 and, in addition, 
the stress was high enough to cause slip in all the grains. The 
experimental distribution is plotted in figure 9 along with the face- 
centered-cubic frozen distribution and the single-mode distribution 
for R = 2.000 taken from reference 1. The comparison is somewhat 
better between the experiment and the face-centered cubic theory than 
between the experiment and the single -mode theory. At the higher slip 
angles, however, the agreement is still not very close. In particular. 
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the maximum experimental slip angle of 58 ° is still far below that pre- 
dicted by the face-centered cubic theory (70° 32 f )* It is a point of 
interest that the analysis of other photomicrographs of similar speci- 
mens strained in the same manner as the specimen from which this 
experimental distribution was obtained has yielded no slip lines with 
slip angles greater than 63°. Thus experience has shown that there is 
an utter lack of slip angles between 63° and the theoretical 70 ° 32 T 
in polycrystals that have been strained in tension until all the grains 
have slipped. It might be argued that the reason for this lack is due 
to the small theoretical probability of a grain exhibiting a slip angle 
between, say, 58 ° and 70 ° 32 1 (on the order of 0 . 05 ), but a simple 
calculation shows that the probability of at least one of 123 grains 
(which is the size of the sample from which the distribution shown in 
fig. 9 was obtained) exhibiting a slip angle between 58 ° and 70 ° 32 1 is 

1 - (1 - 0 . 05) 123 


or O. 998 , almost a certainty. Some other reason must therefore exist 
for the aforementioned lack. This reason is undoubtedly connected with 
either the fact that surface grains do not act like interior ones or 
the violation of one of the assumptions. With regard to this latter 
point, the assumption dealing with equality of microscopic and macro- 
scopic stresses seems to be the most likely suspect, particularly since 
it is generally agreed that the stress state in a polycrystal varies 
from grain to grain. (See, for example, ref. 7.) 


CONCLUDING REMARKS 


The differences between theoretical slip-angle distributions 
derived on the basis of polycrystal aggregates made up of single-slip- 
mode crystals and those made up of face -centered cubic crystals are 
found to be appreciable, the difference increasing as the stress is 
increased. 

The face-centered cubic theory agrees somewhat better with experi- 
ment than the single -mode theory, but the errors are still appreciable. 
These errors are probably due to the inadequacy of the assumption 
involving the homogeneity of stress from grain to grain. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., June 20, 1952 
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APPENDIX 


EQUATIONS OF GREAT CIRCLES BOUNDING THE TRIANGLE OF INTEREST 


In order "to obtain analytical results for W(x) the equations of 
the great circles bounding the triangle shown in figure k must be known. 
The approach that is used in this analysis is to find first the general 
equation of a great circle and then to fit the equation to the points 
through which each bounding circle must pass. 

In terms of the coordinates X,P in figure 3, the direction 
cosines of a given direction (or a given point on the sphere from 
which the stereographic projection is obtained) with respect to the 
x-, y-, and z-axes are 


I XJ ly, l z = sin X cos 0, sin X sin 0, cos X 


(Al) 


A great circle can be defined, in one way, as the locus of points 
on the sphere that are 90 ° away from a fixed point called the pole of 
the great circle. If the location of the pole of a great circle is 
denoted by X 0 ,Po with the direction cosines Z Xq , l yQ , 1 Zq , the 

points on the great circle must satisfy the equation 


l x Z x Q + Wo + lzlz 0 ° 


or 


sin X sin Xq cos(P - 0q) + cos X cos Xq - 0 


or 


0 = 0q + cos“^-(-cot X cot Xq) 


(A2) 


Equation (A2) is the general equation for a great circle. The 
equation has two unknown constants, 0 q and cot Xq. These constants 
can be determined from the coordinates of two points through which the 
great circle is known to pass. The corners of the spherical triangle 
shown in figure 4 furnish the necessary points. 


NACA TN 2777 


23 


The coordinates of the corners can be found from__the angles 
between the corners and the known points (ill) and (Oil). Since the 
Miller indices of the various directions are direction numbers referred 
to a common set of axes, these angles can be easily determined. For 
instance, the cosine of the angle between the (ill) point and the 
(101) point is given by 



and the cosine of the angle between the (Oil) point and the (101) point 
is given by 


( 0 ) 




1 

2 


Knowledge of these two angles allows the determination of the coordinates 
of the point (101). Thus 


(sin X cos (3 ) (0) + (sin X sin 0)(O) + (cos X)(l) 



(sin X cos P)(l) + (sin X sin 0)(O) + (cos X)(0) 


1 

2 


therefore. 


cos X = 


and 
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Similarly, the coordinates of the other two corners can he determined 
to he 


cos X = — 


> for the corner (ill) 


3 = 7 


and 


cos A, 


3 = - ■ 


for the corner (001) 


In order that each hounding great circle pass through the required 
points the arbitrary constants in equation (A2) must he such that the 
resulting equations are: 

For great circle (001,101), 


3 


Tt 

6 


cot X 
cos 

\f2 


(A 3a) 


For great circle (001,111), 


p = cos“l (^2 cot \) - -2- 


(A 3b) 


For great circle (101,111), 



(A3c) 


4B 
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Figure 1.- Slip systems of the face-centered, cubic crystal. 




Figure 2.- Stereographic projection for obtaining the most highly loaded 
slip system for a known orientation. The numbers within each region 
indicate the most highly loaded slip system. 
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Figure 3.- Geometry of the planes, directions, and angles involved 

in the analysis. 


28 


NACA TN 2777 



Figure 4.- Stereographic projection of the curve T = T g superimposed 
on the particular spherical triangle of interest. 
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Figure 5.- Slip-angle distribution for stress ratio R - 1.015. 
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Figure 6.- Slip-angle distribution for stress ratio R = 1.156. 
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Figure 7-- Slip-angle distribution for stress ratio R = 1.837. 
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Figure 8.- Variation of maximum slip angle with stress ratio. 



Figure 9.- Theoretical and experimental slip-angle distribution for 

stress ratio R = 2.000. 
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